Reheating- volume measure in the landscape 
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I recently proposed the "reheating- volume" (RV) prescription as a possible solution to the measure 
problem in "multiverse" cosmology. The goal of this work is to extend the RV measure to scenarios 
involving bubble nucleation, such as the string theory landscape. In the spirit of the RV prescription, 
I propose to calculate the distribution of observable quantities in a landscape that is conditioned in 
probability to nucleate a finite total number of bubbles to the future of an initial bubble. A general 
formula for the relative number of bubbles of different types can be derived. I show that the RV 
measure is well-defined and independent of the choice of the initial bubble type, as long as that 
type supports further bubble nucleation. Applying the RV measure to a generic landscape, I find 
that the abundance of Boltzmann brains is always negligibly small compared with the abundance of 
ordinary observers in the bubbles of the same type. As an illustration, I present explicit results for 
a toy landscape containing four vacuum states and for landscapes with a single high-energy vacuum 
and a large number of low-energy vacua. 



I. INTRODUCTION AND SUMMARY 

In many cosmological scenarios the fundamental the- 
ory does not predict with certainty the values of observ- 
able cosmological parameters, such as the effective cos- 
mological constant and the masses of elementary parti- 
cles. This is the case even for some models of inflation 
driven by a scalar field (see e.g. [E for early work) as 
well as for the "landscape of string theory" [E 0, IE @l j see 
also the "recycling universe" Q models. In these latter 
models, the fundamental theory admits a large number of 
disjoint vacuum states. Transitions between these states 
are possible through bubble nucleation; the interior of 
a bubble appears as an infinite homogeneous open uni- 
verse [H , if one disregards the small probability of bubble 
collisions (see Refs. [E Ei| for analyses of bubble colli- 
sions) . The presently observed universe is situated within 
a bubble (called a "pocket universe") of some type. 

A common feature of these cosmological models is 
the presence of eternal inflation, i.e. the absence of 
a global end to inflation in the entire spacetime (see 
Refs. 0E [l3| for reviews). Eternal inflation gives rise 
to infinitely many causally disconnected regions of the 
spacetime where the cosmological observables may have 
significantly different values. In the context of the string- 
theoretic landscape, eternal inflation entails the nucle- 
ation of (potentially) infinitely many nested bubbles of 
different vacuum types. 

The program outlined in the early works [IE [IE Eil j 
which dealt with eternal inflation of random walk type, 
was to calculate the probability distribution of the cos- 
mological parameters as measured by an observer ran- 
domly located in the spacetime. The main diffuculty in 
obtaining such probability distributions is due to the infi- 
nite volume of regions where an observer may be located. 
An eternally inflating universe contains an infinite, inho- 
mogeneous, and topologically complicated spacelike hy- 
persurface (the reheating surface) where observers may 
be expected to appear with a constant density per unit 
3- volume. 



In the landscape scenarios, one encounters a kind of 
infinity that is in some sense more ill-behaved than in 
the random-walk inflationary scenarios. Not only each 
pocket universe may contain infinitely many observers, 
but also the number of different pocket universes in the 
entire spacetime is infinite. Pocket universes of differ- 
ent types are not statistically equivalent to each other 
because they have different rates of nucleation of other 
pocket universes. There seems to be no natural ordering 
on the set of all pocket universes throughout the space- 
time, since most of the pocket universes are spacelike sep- 
arated. To emphasize the mutual causal independence of 
pocket universes, one calls such a spacetime a "multi- 
verse." 

In summary, eternal inflation is a stochastic process 
that generates a topologically complicated and noncom- 
pact locus of points where observers may appear. A 
"random location" of an observer within that locus is a 
mathematically undefined concept, similarly to the con- 
cept of an integer number "uniformly chosen" among 
all the integers, or a real number "uniformly chosen" 
among all the reals. This is the root cause of the 
technical and conceptual difficulties known collectively 
as the "measure pro blem" in multiverse cosmology (see 
Refs. [H EE El EE EE 13 for reviews). Nevertheless, 
one may try to formulate a prescription for calculating 
probabilities of observer-based events. Such a prescrip- 
tion, also called a "measure," should in some sense corre- 
spond to the intuitive notion of probability of observation 
at a "random" location in the spacetime. 

Several measure prescriptions have been proposed in 
the literature. The proposals that apply directly to land- 
scape scenarios are the "holographic" measure [IE [22[ 
(see also the recent proposal (23[), the "como ying horizon 
cutoff" [H HE SI, the "stationary measure" [IE S3, the 
measure on transitions [28II. and the "pseudo-comoving" 
measure [IE HE HE l32j|. In the absence of a unique 
definition of the measure, one judges a cutoff prescrip- 
tion viable if its predictions are not obviously patho- 
logical. Possible pathologies include the dependence on 
choice of spacetime coordinates [IE HH, the "young- 
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ness paradox" [3 
lem [2i[3l[3i,l 



and the "Boltzmann brain" prob- 
iHEEEIj]- Various measures have 
been used for predicting cosmological parameters, most 
notably the cosmological constant, in the landscape sce- 
narios (see, e.g., Refs. [13, H EE EE S3, SE S3| ) - 

The purpose of this paper is to extend the most re- 
cently proposed "rcheating-volume" (RV) measure (50l . 
5lj . originally formulated in the context of random- walk 
inflation, to landscape scenarios. The basic idea of the 
RV proposal is to select multiverses that are very large 
but (by rare chance) have a finite total number of ob- 
servers. In the context of a string landscape scenario 
(or a "recycling universe"), this can happen if sufficiently 
many anti-de Sitter or Minkowski bubbles nucleate every- 
where, collide, and merge. In such a case, there will be 
a finite time after which no de Sitter regions remain and 
no further nucleations can occur. Hence, there will be a 
finite time after which no more observers are created any- 
where. By this construction, one obtains a subensemble 
of multiverses having a fixed, finite total number 7V b s of 
observers. These finite multiverses with very large iV b s 
are regarded as controlled approximations to the actual 
infinite multiverse. The observer-weighted statistical dis- 
tribution of any quantity within a finite multiverse can 
be obtained by ordinary counting, since the total num- 
ber of observers within any such multiverse is finite. The 
limit of that statistical distribution as N b s — > oo is the 
final result of the RV prescription. 

It was shown in Refs. [HE HH that the RV measure 
is gauge-invariant, independent of the initial conditions, 
and free of the youngness paradox in the context of 
random-walk inflation. Presently I investigate whether 
the same features persist in an application of the RV mea- 
sure to landscape scenarios. In particular, it is important 
to obtain RV predictions with respect to the "Boltzmann 
brain" problem that has been widely discussed. 

In principle, the RV prescription can be extended to 
landscape models in different ways, depending on the 
precise choice of the ensemble of finite multiverses. The 
ensemble of multiverses with a fixed total number of ob- 
servers N bs (where one counts both the ordinary ob- 
servers and the "Boltzmann brains") appears to be the 
natural choice. However, it is difficult to compute the 
number of observers directly and unambiguously. Instead 
of the total number of observers, I propose to fix the to- 
tal number n to t of bubbles nucleated to the future of an 
initial bubble. 

The total number of ordinary observers in bubbles of a 
given kind is proportional to the volume of the reheating 
surface in those bubbles. It is known that the "square 
bubble" approximation [42|, which neglects the effects of 
bubble wall geometry, is adequate for the purposes of vol- 
ume counting. Then the evolution of the landscape is well 
described by the approximate model called "inflation in a 
box" [HE H2, [SU • In that approximation, one keeps track 
only of the number of new bubbles nucleated in previ- 
ously existing bubbles, and each new bubble is assumed 
to be instantaneously nucleated exactly of Hubble size in 



comoving coordinates. Motivated by this approximation, 
in this paper I study a simplified definition of the RV mea- 
sure for a landscape scenario (see Sec.|II]for details): One 
requires the total number of bubbles of all types, n to t, to 
be finite and evaluates the statistical distribution of bub- 
ble types (or other cosmological observables) in the limit 
ntot ~~ * °o- I n principle, this limit can be calculated if the 
bubble nucleation rates are known. Since this paper is a 
first attempt to perform this techically challenging cal- 
culation, I concentrate only on bubble abundances and 
on the relative abundance of Boltzmann brains. I neglect 
the increased number of observers due to additional slow- 
roll inflation within bubbles; this effect was considered in 
Ref. [!(| and requires additional complications in the for- 
malism. 

A landscape scenario may be specified by enumerating 
the available N types of vacua by a label j (j = 1, N) 
and by giving the Hubble rates Hj within bubbles of type 
j. One can, in principle, compute the nucleation rate 
Tj^k describing the probability (per unit four-volume of 
spacetime) of creating a bubble of type k within bubbles 
of type j. It is convenient to work with the dimensionless 
rates, 
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(1) 



The rate Kj_>fc equals the probability of having a bubble 
of type k within a 3- volume of one horizon in a bubble of 
type j, during a single Hubble time. 2 Explicit expressions 
for available in some landscape scenarios. In 

what follows, I assume that known. 

In Sec. IIIII I apply the RV measure proposal to a toy 
model of the landscape with four vacua (the FABI model 
of Ref. [13] )■ I n this model, the vacua labeled F and I 
are de Sitter (dS) and the vacua labeled A and B are 
anti-de Sitter (AdS) states. One assumes that only the 
transitions F —> I, I —> F, F —> A, and I — > B are 
allowed, with known nucleation rates Kfi, Kif, etc., per 
unit Hubble 4-volume. I show in Eq. (fr?3")) that the RV- 
regulated bubble abundances depend on the value of the 
dimensionless number 



n 



\ u+l 

kib \ kfi 



v = e 3 . 



UFA, 



KiF 



(2) 



Here the constant v, introduced for convenience, is simply 
the number of statistically independent Hubble regions 
after one e-folding. Barring fine-tuned cases, one expects 
that the value of rj is either much larger than 1 or much 
smaller than 1, since the nucleation rates may differ by 
exponentially many orders of magnitude. By relabeling 
the vacua [F <-> I and A <-> B) if necessary, we may 



1 For some recent work concerning the determination of the bubble 
nucleation rates, see Refs. |54l. l55l l56l |57|. |5£| . 

2 The notation K,_>fc, chosen here for its visual clarity, corresponds 
to Kfej of Ref. |24| and to r,-(, of Ref. [Hol l. 
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assume that 77 <C e 6 . Then the bubble abundances are 
approximately described by the ratios 



p(I) : p(F) : p(A) : p(S) » \ : i : 1 



n 



(3) 



This result can be interpreted as follows. Each of the I 
bubbles produces v bubbles of type F, and each of the F 
bubbles produces v bubbles of type A. The abundance 
of B bubbles is neligible compared with other bubbles. 
Heuristically, the chain of transitions I — > F — > A can 
be interpreted as the "dominant" chain in the landscape. 
The fine-tuned case, e -6 < r\ < e 6 , is considered sepa- 



rately, and the result is given by Eq. (|54j) . 

I then consider the abundance of Boltzmann brains 
(BBs) in the FABI landscape (Sec. iHTBj) . The total 
number of BBs is proportional to the total number of 
Hubble regions (if-regions, or 4-volumes of order i7~ 4 ) 
in de Sitter bubbles. The coefficient of proportionality is 
the tiny nucleation rate r BB of Boltzmann brains, which 
is of order exp(— 10 50 ) or smaller. In comparison, ordi- 
nary observers occur at a rate of at least 1 per horizon 
volume. It turns out that (after applying the RV cut- 
off) the total number of £f-regions of dS types F or I 
is approximately equal to the total number of nucleated 
bubbles of the same type. Hence, the BBs are extremely 
rare compared with ordinary observers. 

In Sec. IIVI I extend the same calculations to a general 
landscape with an arbitrary number of vacua. The RV 
prescription predicts a definite ratio p(j)/p(k) between 
the number of bubbles of types j and k. I derive a for- 
mula for the ratio p(j)/p(k) that involves all the param- 
eters of the landscape. With the help of mathematical 
results derived in Sections IIV El and IIVFI it is possible 
to show in full generality that the RV measure gives well- 
defined results that are independent of the initial condi- 
tions. Nevertheless, actually performing the required cal- 
culations for an arbitrary landscape remains a daunting 
task. To obtain explicit expressions in a semi-realistic 
setup, I calculate the bubble abundances for a landscape 
that contains a single high-energy vacuum and a large 
number of low-energy vacua. The result is an approxi- 
mate formula [Eq. (|124|l ] for the ratio p(j)/p{k) expressed 
directly through the nucleation rates of the landscape. 

I also demonstrate in Sec. IIV CI that the abundance of 
Boltzmann brains is negligible compared with the abun- 
dance of ordinary observers in the same bubble type. 

To conclude, the present work demonstrates that the 
RV measure has attractive features and may be consid- 
ered a viable candidate for the solution of the measure 
problem in multiverse cosmology. More work is needed to 
investigate the dependence of the predictions on the pre- 
cise details of the definition of the ensemble E n . I have 
developed an extensive mathematical framework for the 
calculations in the RV prescription and obtained first re- 
sults for specific landscapes. However, a more powerful 
approximation scheme is desirable so that the predictions 
of the RV measure can be more easily obtained for land- 
scapes of general type. Ultimately, the viability of the 



RV measure is to be judged by its predictions for cosmo- 
logical observables in realistic landscapes. These issues 
will be considered in future publications. 



II. THE RV MEASURE FOR A LANDSCAPE 

The RV measure prescription as formulated in 
Refs. HU applies only to the calculation of abun- 
dances of terminal bubbles. We will now extend the RV 
measure to computing arbitrary statistics on a landscape. 

We first note that RV measure prescription can be ap- 
plied, strictly speaking, only to landscapes that contain 
some terminal bubble types (i.e. vacua from which no 
further tunneling is possible). However, this limitation is 
quite benign, for two reasons. First, a landscape with- 
out any Minkowski or AdS states is not expected to be 
realized in any realistic string theory scenario without 
an exceptional amount of fine-tuning. Second, the pre- 
viously proposed volume-based and the world-line based 
measure prescriptions agree for a landscape without ter- 
minal bubble types (2ll |26|. One may therefore consider 
the measure problem as solved in such landscapes and 
turn one's attention to more realistic landscapes where 
terminal bubble types are present. 

Let us take an initial bubble of a nonterminal type j 
and consider the statistical ensemble E n (j) of all possible 
evolutions of the initial bubble such that the total number 
of nucleated bubbles of all types is finite and equals n 
(not counting the initial bubble). The total number of 
nucleated bubbles in a multiverse can be finite only if 
terminal bubbles nucleate everywhere and merge globally 
to the future of the initial bubble. This can happen by 
rare chance; however, it is important the total probability 
of all events in the ensemble E n (J) is always nonzero for 
any given n, so that the ensembles E n (j) are well-defined 
and nonempty. 

The ensemble E n (j) may be described in the language 
of "transition trees" used in Ref. [2l[ . The ensemble con- 
sists of all trees that have a total number n+ 1 of bubbles, 
including the initial bubble of type j. The trees in E n (j) 
are finite because all the "outer" leaves are bubbles of ter- 
minal types. The motivation for considering the ensemble 
E n (J) is that a finite but very large tree (with n large) is 
a controlled approximation to infinite trees that typically 
occur. Hence, we are motivated to consider E n (j) with 
n finite but very large. 

Note that the ensemble E n (j) differs from the ensemble 
defined in Ref. [50]; in E n (j) the total number of bubbles 
of all types is equal to n, rather than the total number 
of terminal bubbles as in Ref. 50]. Thus, the current 
proposal, which appears more natural, is an extension 
of that of Ref. [50( . Future work will show whether this 
technical difference is significant; presently I will investi- 
gate the consequences of the current proposal. 

Once the ensemble E n (j) is defined, one may consider 
the statistical distribution of some cosmological observ- 
able within the multiverses belonging to the set E n (j). 
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For instance, one can count the number of bubbles of 
some type k, or the number of observers within bubbles 
of type k, or the number of observations of some physi- 
cal process, etc. In a very large multiverse belonging to 
E n (j) with n 3> 1, one may expect that the statistics 
of observations will be independent of the initial bubble 
type j. Indeed, this will be one of the results of this 
paper. Hence, let us suppress the argument j and write 
simply E n . 

Each multiverse belonging to E n has a naturally de- 
fined probability weight, which is simply equal to the 
probability of realizing that multiverse. This probability 
weight needs to be taken into account when computing 
the statistical distributions of observables. The sum of all 
probability weights of multiverses within E n is equal to 
the total probability of E n , which is exponentially small 
for large n but always nonzero. Since the multiverses 
from the ensemble E n are by construction finite, i.e. each 
multiverse supports only a finite total number of possible 
observers, we are assured that any statistics we desire to 
compute on E n will be well-defined. 

We can now consider the probability distribution 
p(Q\E n ) of some interesting observable Q within the en- 
semble E n and take the limit n — > oo. One expects that 
the probability distribution p(Q\E n ) will have a well- 
defined limit for large n, 

p(Q) = lim p{Q\E n ). (4) 

n — >oo 

It was shown in previous work on the RV prescrip- 
tion [HclHH that the distribution p(Q) is well-defined for 
a simplest toy landscape as well as in the case of random- 
walk eternal inflation. In this work I extend these results 
to a general landscape scenario. Below (Sec. HVEp I will 
prove rigorously that the limit ([J| indeed exists and is 
independent of the chosen initial bubble type j as long 
as the initial bubble is not of terminal type and as long as 
the landscape is irreducible (every vacuum can be reached 
from every other non-terminal vacuum by a chain of nu- 
cleations). Thus, the distribution p(Q) is unique and 
well-defined. This distribution is the final result of ap- 
plying the RV prescription to the observable Q. 

In practice, it is necessary to compute the distribu- 
tion p(Q\E n ) asymptotically in the limit of large n. A 
direct numerical calculation of probabilities in that limit 
by enumerating all possible evolution trees is extremely 
difficult because of the exponential growth of the number 
of possible evolutions. Instead, I derive explicit formulas 
for the distribution p(Q) in a generic landscape by eval- 
uating the limit (0| analytically. These formulas are the 
main result of the present article. 

III. A TOY LANDSCAPE 

I begin by applying the RV measure to a toy landscape 
with very few vacua. Using this simple example, I de- 
velop the computational techniques needed for the prac- 
tical evaluation of the limit such as Eq. In Sec. IIVI 



the same techniques will be extended to a more general 
landscape with an arbitrary number of vacua. 

In Ref. I considered the simplest possible nontrivial 
landscape: a single dS and two AdS (terminal) vacua. 
The next least complicated example that can be treated 
analytically is a toy landscape having two nonterminal 
and two terminal vacua. This toy landscape was called 
the "FABI" model in Ref. [24| and consists of the vacua 
labeled F, I, A, B with the transition diagram A <— F «-> 
I — > B. In other words, one assumes that the vacuum F 
("false" vacuum) can nucleate only bubbles of types A and 
i, the vacuum I ("intermediate" vacuum) can nucleate 
bubbles of types F and B, while A and B are terminal 
vacua that do not have further nucleations. 

To describe the finitely produced probability in this 
landscape, I use the discrete picture called the "eternal 
inflation in a box" [Hj], which is closely related to 
the "square bubble" approximation [42| . In this picture, 
one considers the evolution of discrete, causally disjoint 
homogeneous iJ-regions in discrete time. All possible 
vacuum types are labeled by j = 1,...,N. During one 
time step of order St — HJ 1 , where Hj is the local Hubble 
rate in a given ii-region of type j, the evolution consists 
of expanding the ii-region into 

e 3fl>« = e 3 _ v (5) 

daughter ii-regions of type j. Each of the daughter H- 
regions has then the probability nj^k of changing imme- 
diately into an if-region of type k =/= j; this imitates a 
nucleation of a horizon-size bubble of type k. If no tran- 
sition has taken place, the daughter ii-region retains its 
type j. For convenience, we denote by 

Kj^j = 1 - K 3^k (6) 
Mj- 
the probability of no transitions during one Hubble time. 
The process of expansion and transition is continued 
ad infinitum, independently for each resulting ii-region. 
The ii-regions of terminal types will admit no further 
transitions and will not expand further (except, perhaps, 
by a fixed amount due to slow-roll inflation occuring im- 
mediately after nucleation). An example simulation is 
shown in Fig. [T] 

In the remainder of this section I perform explicit cal- 
culations of the RV cutoff in the FABI model. 



A. Bubble abundances 

The first task is to compute the relative abundance 
of bubbles of different types. Consider the probability 
p(n to t, n>F,ni, ua, ns',k) of having a finite total number 
ntot of bubbles of which nj are of type j (where j = 
F, I , A, B), if one starts from a single initial ii-region 
of type k (k = F,I). By construction, p ^ only for 
That — n F + Jfy + tia + n-B- A generating function for this 
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Figure 1: An example simulation of "eternal inflation in a 
box" in two spatial dimensions. Bubbles (_ff-regions) are rep- 
resented in comoving coordinates by squares. Dark shades 
indicate bubbles of terminal types. Other shades correspond 
to nested bubbles of various nonterminal ("recyclable") types. 
For the purposes of visual illustration, nucleation rates were 
chosen of order one, and colored lines were drawn at bubble 
boundaries. 



probability distribution can be defined by 

g(z,{qj};k) : 



£ z>(n,K};fc) JJ 
n>0,n 3 >0 j=F,I,A,B 

(z n &tfqr<&°) n<0B . ik , (7) 



where the notation (— } n<0O .^ stands for a statistical av- 
erage restricted to events with a finite total number n 
of bubbles nucleated to the future of an initial bubble of 
type fc. 

The generating function g plays a crucial role in the 
entire calculation. Since we will be only interested in the 
initial bubbles of types F and I, let us denote 

F(z, {qj}) = g(z, {<?,} ; F), I(z, { Qj }) = g{z, {q,} ; 7). 

(8) 

The generating functions F and I satisfy the following 
system of nonlinear algebraic equations [50j , 

F« = zq A n FA + zqiKFiI + n FF F, (9) 
I~ = zq B n IB + zq F n IF F + Kul. (10) 

Here we denoted for brevity k ff = 1 — K F a — kfi and 
kjj = 1 — kjb — kj F ; within our assumptions, k ff w 1 
and kij w 1. In Eqs. (f9"))- (fTU|) the generating variable 
z multiplies only the terms that correspond to changing 



the type of the iJ-region (which imitates the nucleation 
of new bubbles) but not the terms k ff F and k.hI that 
correspond to the eventuality of not changing the type of 
the iJ-region during one Hubble time. 

The nonlinear equations (|9|)— (p~0|) may have several 
real-valued solutions, as well as complex- valued solutions 
that are certainly not of physical interest. In particular, 
for z = 1 and qj = 1 there exists the "trivial" solution 
F = I = 1 as well as a nontrivial solution with F <C 1 and 
J « 1. Similarly, for z near there exists the solution 
that approaches F(0) = 1(0) = 0, 

F = (zq A K FA y + Oiz 2 - 1 ), I « {zq B K IB y + Oiz 2 "- 1 ), 

(11) 

as well as the solution 



F : 



"FF 



1, J: 



"II 



1 



(12) 



and solutions where F«l and / sa and vice versa. It 
is important to determine the solution branch F{z), I(z) 
that has the physical significance as the actual generat- 
ing function of the finitely produced distribution of H- 
regions. 

The functions F(z),I(z) are solutions of algebraic 
equations and thus are continuous functions of z that 
are analytic everywhere in complex z plane except for 
branch cuts. It is easy to see that the solution F = I = 1 
at z = 1, qj = 1 is continuously connected with the solu- 
tion (fT2"|) at z ~ 0, while the solution (fTl"]) is continued to 
a solution with F(z) <C 1 and I(z) <C 1 for all < z < 1. 
The values F(l) and 1(1) are the probabilities of the 
events that the evolution of an initial bubble of type F or 
I ends globally. These probabilities are extremely small 
and of order k fa and k ib respectively. This is easy to 
interpret because, for instance, k fa is the probability of 
nucleating v terminal regions at once after one Hubble 
time within an 77-region of type F. Hence, the gener- 
ating functions for the finitely produced distribution of 
bubbles are given by the solution branch having F <C 1 
and J « 1 rather than by the solution F = I = 1 at 
z = 1. More precisely, the physically meaningful solution 
F(z),I(z) is selected by the asymptotic behavior (fTl"]) . 
This argument removes the ambiguity inherent in solving 
Eqs. (|5|)- (fTU)) . Below we refer to the branch of solutions 
F(z),I(z) connected to the nontrivial solution pT|) near 
z = as the "main branch." 

Once the main branch of the generating functions 
g(z 1 {qj}] k) are known, one can express the mean num- 
ber of bubbles of type i (i = F, I, A, B) at a fixed total 
number ritot as 



p{i\n to t) 



d^d qi g(z = 0,{ qj = l};k) 



That n tot dT at g{z = 0, {qj = 1} ; fc) ' 

One expects that the limit of this ratio at ntot - * oo 
will be independent of the initial bubble type fc since the 
ensemble E ntat will consist of _ff-regions having a very 
long evolution, so that the initial conditions are forgotten. 
Below I will show explicitly that this is indeed the case. 
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It is more convenient to compute the ratios of the num- 
ber of bubbles of types i and i' at fixed n tot , 



p(z>tot) d^d^g 



(14) 



z =0,{ qj = l} 



Then one only needs to compute derivatives d qi g = g. qi 
evaluated at qj = 1. These derivatives satisfy a system of 
linear equations that can be easily derived from Eqs. 
(fTU|) . For instance, the derivatives F <qA and I iQA satisfy 



1 



- 1 / 



9-1 



ZK F A + ZK FI I. qA + K FF F tqA , 
ZK IF F^ qA + 

IS 

Rewriting these equations in a matrix form, we obtain 

±-Fi^-K FF 



(15) 
(16) 



-ZKjp 



-1^ 



ZKFI \ 


c ,QA 




ZKFA 


1 - Kll ) 


ilA 








The coefficients of the ^-dependent matrix 

^Fi^ 1 - K F F -ZKpi 
1 r^-l 

-ZK IF -I" - Ku 



M(z) 



(17) 



(18) 



are the main branch of solutions of the nonlinear equa- 
tions (j9"))- (fl"t)|) at qj — 1 but at arbitrary z. 

One can verify using Eq. (jTT|) that the matrix M(z) is 
invertible near z — 0. Hence the solution of Eq. (fT7|) can 
be written, at least within some range of z where M(z) 
remains invertible, as 



■'I A 

1A 



M-\z) 



ZKFA 





(19) 



Similarly, the derivatives Fq F and I iqF satisfy the equa- 
tions 



M{z) 

whose solution is 
p 










>1 p 




ZKjpF 



= M-\z) 





zkifF 



(20) 



(21) 



Other generating functions can be expressed in the same 
manner. 

One could in principle obtain a numerical solution for 
F tqA {z), I tqA {z), and all the other generating functions 
at any given value of z. However, the numerical solution 
is not particularly useful at this point because the next 
step in the calculation is the evaluation of Eq. (fT4"|) in the 
limit of very large ntot, for instance, 



p(A\n tot ) _ 8?"F t 
p(F\n tot ) 



QA 



(22) 



It is generally not feasible to compute the n-th derivative 
of a numerically obtained function in the limit n — > oo, 



because the unavoidable round-off errors are amplified by 
a fixed factor with each successive numerical differentia- 
tion. Therefore, we need a way to evaluate the ratio ([T4]) 
in the limit n t ot ™ ► oo without using numerics. Indeed 
we will be able to compute the ratio p(A)/p(F). We will 
also show that, for instance, 



P(A) = Hm d™F qA 
p(F) n—»x 



d?F 



IF 



f) n T 
= lim ° z ' qA 



z=0 



d?i. 



'IF 



(23) 



z=0 



in other words, that the final RV-regulated ratio 
p(A)/p(F) is independent of whether the initial bubble 
is of type F or of type /. 

To proceed, we use the fact that F tqA (z), I 9A (z), 
F^ qF (z), etc. are analytic functions of the parameter z. 
The asymptotic growth of high-order derivatives of an 
analytic function f(z) is determined by the location of 
the singularities of f(z) in the complex z plane. The re- 
quired result can be derived by the following elementary 
argument: Consider the derivative d n f/dz n at z = 0, 
and assume that f(z) admits an expansion around the 
singularity z* nearest to z — 0, such as 



f(z) = c + ci (z - z„) s + 



(24) 



where s 0,1,2,... is the power of the leading-order 
singularity, and the omitted terms are either higher pow- 
ers of z — z* or singularities at points z'^ located further 
away from z — 0. The singularity structure (|24[) yields 
the large-n asymptotics 



dz n 



,r(n 



~c 1 (-z») s 
z=o r (- s ) 



(25) 



It can be seen from this formula that any other singular 
point zl located further away from z = gives a con- 
tribution that is smaller by the factor |z^/z*| s ™. The 
contribution of a subdominant singularity of the form 

(z — z*) s , i.e. at the same point z = z„ but with a higher 
power s' > s, is suppressed, in comparison with the term 
in Eq. (|23|) . by the factor 



r(- 



T(n-s') r(-s) _ 1 
r(n-a) r(-a') ^ („ _ if- s r(-s') ' 



(26) 



It is clear that the terms omitted from Eq. (|2"5"]) indeed 
give subleading contributions at large n, and so Eq. (j2"5|) 
is indeed the leading asymptotic term at n — * oo. 

We now need to determine the location of the singulari- 
ties of F tQA (z), I, qA (z), etc., as functions of z. This task is 
much simplified once we observe that all these quantities 
are expressed through the inverse matrix M _1 (z), and 
hence it remains to analyze the singularities of M _1 (z). 
That matrix can be singular at some value z — z* cither 
because some of the coefficients of M(z) are singular, or 
because the matrix M(z) is degenerate (noninvertible) 
at z — z+. The coefficients of M(z) depend on solutions 
F(z),I(z) of algebraic equations (fU|)— (fTUj) at qj — 1 and 
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thus cannot be divergent as functions of z. These coeffi- 
cients can be singular only in that some derivative in z 
diverges. The derivatives F 2 and I. z satisfy the equations 



M(z) 







Ufa ~ 


h K FI I 






KlB ^ 


- kifF 



(27) 



Therefore, F x and I z diverge only for those z for which 
the matrix M(z) is degenerate. For these z, all the deriva- 
tives F t g, I tZ , F t q A , etc. will be divergent at the same 
time since they are all proportional to the inverse matrix 
Af-^z). 

We note that the matrix M(z) is the Jacobian of the 
nonlinear system 1[9])-(|1Q|), As long as M[z) is nondegen- 
erate, all the different branches of the solutions F(z), I(z) 
do not meet and remain smooth functions of z. As we 
noted above, the "main" branch F(z), I(z) is the one con- 
nected to the nontrivial solution (fTTj) at z w 0. The ma- 
trix M(z) is nondegenerate near z — 0; therefore, the 
solutions F(z),I(z) remain smooth functions of z for all 
z such that M(z) is nondenegerate. We conclude that the 
only possible singularities of F(z),I(z) are those values 
z = z* where det M(z) = 0. 

Below it will be shown that the behavior of detM(z) 
near z = z* is 



det M(z) w ciV^* — z , 



(28) 



i.e. of the form (f2"4"|) with cq = and s — \. It now 
follows from Eqs. (fl9|) and (|21jl that the derivatives such 
as F >QA (z), I iqA (z), F tqp (z), etc., all diverge at z = z* 
with the same asymptotic behavior, namely proportional 



to (z* — z) 



-1/2 



M-!(z) = 
M(z) = 



. We may express the inverse matrix as 
' (29) 

, (30) 



-M(z), 
det M(z) 

1 7-^-1 

-J" - K U ZK F I 

zk if i-Fi- 1 - 



K F F 



where we introduced the algebraic cofactor matrix M = 
M~ det-M, which is nonsingular at z = z*. It then 
follows from Eq. (fTi?)) that the asymptotic behavior of 
the functions F qA (z) and J qA (z) near z = z* is given by 



F 



1A 



ci^/zl 



=M(z* 



Z*K F A 





C\\JZ* 



-I» l -kh)z*kfa 
zIkif^fa 



(31) 



Similarly, using Eq. (|2"Tj) we find near z = z» 



■</F 



1 



1 



M(*„) 





z^kfikifF 

(r;^" -1 _ K ^ 1 z * k ifF 



(32) 



Using Eq. 



and noticing that all n-dependent factors 



are the same in d z F^ qA and other such derivatives, we 
conclude that 

U~ r .a 



lim 



B n F 



i- 



KFA 



lim ^ ' 9A 



(9"/ 



z^kfiKifF* 
z*k F a 



- K F F F» 



(33) 
(34) 



It is important that the two limits above are equal; this 
is so because the condition det M{z*) = yields 



detM 



k F f 



it 1 



-z^kfik-if 



0. 



(35) 

where we denoted F* = F(z*), i* = I(z*) for brevity. 
It follows that the ratio of yl-bubbles to F-bubbles is 
independent of whether the initial bubble is of type F 
or of type /. In other words, the RV-regulated ratio of 
^4-bubbles to F-bubbles is 



P(A) 
p(F) 



z*kfa 



kff F* 



(36) 



independently of the initial bubble type. Below 
(Sec. IIVE[) the independence of initial conditions will be 
rigorously proved for a general landscape using mathe- 
matical techniques of the theory of nonnegative matri- 
ces. Presently we have shown this independence using 
explicit formulas available for the toy landscape under 
consideration. 

In a similar way, we find the RV-regulated ratio of A- 
bubblcs to 5-bubbles, 



p(A) 



z*kfak-if 



^ (±Fr l KFF) K IB 

and the ratio of F-bubbles to /-bubbles, 
p(F) _ (h< 



(37) 



■-F2 



KFF ) F t 



P(I) 



ZtfKFlI* 



(38) 



It remains to compute z*, F*, i* and to justify Eq. (|28|) . 
We will do this using the explicit form of the matrix 
M(z). To determine z*,F*, and /*, we need to solve 
Eq. (|33|) simultaneously with the equations 

= z*k F a + z*kfiI* + kffF*, (39) 

i_ 

I* = Z*KIB + Z«KifF« + Kill*: (40) 

the latter two being Eqs. (f9"|)- (fTU)) after setting qj = 1. 
We are interested in the solution z* closest to z = 0. By 
definition l[7|). the generating functions are nonsingular 
for |z| < 1; hence, the only possible values of z* are in 
the domain \z\ > 1 in the complex plane. 
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Moreover, we can show that the solutions F(z),I(z) 
are growing functions of z for real z < z*. Initially at 
z w these functions have the form (fTTj) and hence are 
positive. The determinant det M shown in Eq. (j3"5)) is 
also positive for z < z* . One can then find from Eq. (|50|) 

that the algebraic cofactor matrix M and hence the in- 
verse matrix M _1 has all positive elements for those z. It 
is then evident from Eq. ((2ZJ) that d z F > and dj > 
as long as the values of F and I are themselves posi- 
tive. Therefore, the solutions F(z),I(z) are positive and 
growing functions of z as long as detM(z) > 0. 

We may thus visualize the behavior of det M(z) as z 
grows. Both F{z) and I(z) will grow with z, so that the 
terms in square brackets diminish while the second term, 
z 2 kfikif, grows. Eventually the product of the square 
brackets in Eq. (|3"5")l will be balanced by the second term, 
and the determinant will vanish. We need to determine 
the smallest value z — z* for which detM(z) = 0. 

Since the physically significant branch of the solution 
involves always very small values F{z) and I{z) for |z| < 
1, while other (unphysical) branches have either F or 
/ approximately equal to 1, it is reasonable to assume 
that F(z„) <§; 1 and /(z*) <C 1 also at z = z* (the self- 
consistency of this assumption will be confirmed by later 
calculations). Then the terms kffF* ps i* 1 , and ~ 
2» can be disregarded in comparison with fV v and 
in Eqs. (f3"9")) - ([4*0|) . In this approximation, we can simplify 
Eqs. (|39 ]) -([4"0] ) to 



F* = z^kpa + z^kfiI*, 

Ii' = Z^KJB + Z*K IF F*, 



(41) 
(42) 



while Eq. (|35[) becomes 

(kfa + kfi!*) {k-ib + k/f-F 1 *) = v 2 k fi ki F FJ*. (43) 



The ratios (j3"6")) - (f3"5)) are also simplified and can be writ- 
ten more concisely as 

p(A) : p(I) : p(F) : p(B) = k FA : {k F iQ 

Ufa + kfiI* kfa + K F ih 
■ : Kib p ■ (44) 

To determine z» , we will obtain an explicit approximation 
for the main branch F(z), I(z) for all z. For small enough 
z, the solutions of Eqs. (j39|) - (|40)) are well approximated 
by Eq. (HU), 



F(z) fa (zk f a) u , I{z) ~ {zkibT ■ 



(45) 



These solutions are obtained under the assumption that 
the terms zkfa and zkib are numerically small (zkfa *C 
1, zkib < 1) and yet dominant in Eqs. ([39|)-(|4U]). These 
terms only remain dominant as long as 



ufa > kfiI(z), k/s > kifF(z) 



(46) 



Substituting Eq. (TID for F(z) and I(z) into Eq. flUD, we 
obtain the conditions 



z <; k 



_l I KFA 
V K FI 



IB 



Z < 



-l / kib 
k fa I 



\Kif 



(47) 



We need to check whether det M(z) could vanish already 
for some z* within the range (f47|) . Using Eq. (|43|) in the 
regime (|46|) . we find 



[kfaKibY 1 v 2 k F ik,if 



(48) 



This value is within the range (|47|) only if the following 
simultaneous inequalities hold, 



1 



KJB \ 



kfa 



^ « v 2 = e e . 

KIF 



Let us denote by 77 the quantity in Eq. |49|) . 



V = 



\ v+l 

kib \ kfi 



kfa . 



KIF 



(49) 



(50) 



then Eq. (|4*9"|) becomes simply |ln?y| < 6. One would ex- 
pect that 77 is generically either very large or very small, 
so the inequalities (|4"9"]l can hold only in a fine-tuned land- 
scape. Additionally, we need to require 



z* <g; min 



which entails 



kib 

KIF 



KFA Kib 



kfa v_ 
kfi ^/rj' 



(51) 



(52) 



However, this requirement is weaker than Eq. (|49[) since 
typically kfa <C kfi and kib <^ kif, so 



Z* K 



IB 



1 I kfa 
kfi 



-1 -1 1 " _i 



FA- 
(53) 

Let us assume, for the moment, that Eqs. (|49|) and (|52|) 
hold. Then we use Eq. (f44|) to obtain 



p(A) : p(J) : p(F) : p(B) « 1 



07 . 1 



0j. (54) 



Due to the fine-tuning assumption (|49|) . these ratios are 
all within the interval |V~ 2 , v 2 \ = [e~ 6 , e 6 ] . 

Having considered the fine-tuned case, let us now turn 
to the generic case. Generically one would expect that 
the quantity 77 is either extremely large or extremely 
small, and in any case outside the logarithmically narrow 
range [e~ 6 ,e 6 ]. In that case, one of the inequalities ([49]) 
does not hold; generically, either 77 <C v~ 2 or 77 ^> v 2 . It 
follows that det M(z) ^ for all z within the range (|4"T)) . 
To reach the value z* at which det M(z*) = 0, we need to 
increase z further, until one of the terms zkfiI or zkifF 
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in Eqs. (j3"9")) - P0")) becomes dominant and the solution (14"5)) 
becomes invalid. 

It is impossible that both the terms zkfiI and zkifF 
are dominant in Eqs. (|3"9")) - f4T))) at z — z* because then 
Eq. ([43]) would yield a contradiction, 



(55) 



Hence, the determinant detM(z) first vanishes at the 
value z = Zn such that only one of those terms is domi- 
nant in its respective equation. Without loss of general- 
ity, we may relabel the vacua (F <-* /, A <-> B) such that 
the term zkifF becomes dominant in Eq. (|40j) while the 
term zkfa is still dominant in Eq. (|39p . This is equiva- 
lent to assuming r\ <C v~ 2 . In the range of z for which 
this is the case, 

zkfiI(z) <§; zkfa <C 1; zkjb <C zkjfF(z) <§; 1, (56) 

the approximate solution of Eqs. (|3"9" |) - (f4T))) can be written 
as 

F(z) « (z^f ; J(z) « (zkjfF) 1 ' = K V p A K v IF z^ +l > . 

(57) 

With these values of F(z) and /(z), the consistency re- 
quirement (|56[) yields the following range of z, 



< z" < 
This range is nonempty if 



K FA K FI K IF 



V 



KlB 
KFA 



u+l 



KFI 
KlF 



« 1, 



(58) 



(59) 



which is indeed one of the two possible ways that the 
fine-tuning (|49|) can fail. Having assumed that the above 
condition holds, we need to determine the value of z* and 
check that it belongs to the range (|56|) . Using Eq. (|43]) 
in the regime (|56j) . we obtain 



2 .v w — 1 

V Kjf^FI^fA 



(60) 



Substituting this value into the inequalities (|58[) . we find 
the condition 



r? -C -5- < 1, 



(61) 



which holds identically under the current assumption, 



(The relabeling F <-* /, A «-* B is neces- 



sary if the opposite case, 77 » u , holds.) Therefore, z» 



is within the regime (|56p . and our approximations are 
self-consistent, yielding 



F» w 



[^ 2 t/i?KF/] 



KFA 

v 2 kfi 



(62) 



The ratios (|44|) become 
: p(J) : p(F) : 



(63) 



This is the result of applying the RV prescription to a 
generic FABI toy landscape in the second regime. 

It remains to justify the statement of Eq. (|2"8")l . Below 
in Sec. IIVEI I will demonstrate that the property (|28|) 
holds for a general landscape. Here only a simple argu- 
ment is presented to illustrate this property for Eqs. (|4"Tj) - 
(|4"2")l . Using those equations, we can express i* through 
F* and derive a closed algebraic equation for F» , 

F* = z v [k F a + k,fiz u (k ib + kifF*) v ] u = f(z;F»). 

(64) 

The solution F(z) is given by the intersection of the line 
y = F and the curve y = f(z; F) in the y — F plane. The 
function f(z; F) is convex in F for F > 0, z > 0; hence, 
there will be a value z = z* for which the curve /(z*; F) 
is tangent to the line y = F, i.e. /,f(z*; F*) = 1. This 
value of z* will then implicitly determine F*. For values 
F ~ F*, z w z* the dependence of F(z) on z will exhibit 
the singularity behavior of the type a/z* — z. To see this 
formally, we may expand 

F = f(z; F) « F* + f tF (F - F.) + ~/,ff (F - F*) 2 
+ f, z (z - z») . (65) 

Since /j?(z*;F*) = 1, we obtain 



F « F* + 



2/, 2 



FF 



vzT— Z- 



(66) 



This shows explicitly the singularity structure of the 
form (|24[) . The determinant of the matrix M is a smooth 
function of F(z) and I(z) near z = z*. Expressing det M 
as a function only of F, we obtain for z « z* the required 
formula 



det M(z) w (F - F* 



det M tx V-z* - i- (67) 



F» 



B. "Boltzmann brains" 

Let us now use the same techniques to compute the 
relative abundance of "Boltzmann brain" observers to or- 
dinary observers. 

We need to introduce an appropriate set of generat- 
ing functions. "Boltzmann brains" can be created with a 
fixed probability per unit 4-volume, unlike ordinary ob- 
servers who can appear only within a narrow interval of 
time after creation of a given bubble. Let us therefore 
compare the total number of bubbles rij with the total 
number of 4- volumes H~ in bubbles of type j. 

To be specific, let us fix j = F (no Boltzmann 
brains can be expected in a terminal vacuum) . Consider 
the probability p(n to t , if, Np] k) of having a finite total 
number n to t of bubbles of which rip are of type F and 
Np is the total number of ^-regions of type F, if one 
starts from a single initial H -region of type k (k = F, J). 
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The number of "Boltzmann brains" in bubbles of type 
F is proportional to Np with a small proportionality 
constant, Kp B , which we will include at the end of the 
calculation. A generating function g for the probability 
distribution p can be defined by 

g(z,q,r;k)= £ z n q np r Np p(n,n F ,Np;k) 

n,n F ,N F >0 

= (z n q nF r NF ) , . (68) 

For the initial bubble types k = F and k = I, let us 
denote for brevity 

F(z,q,r) = g(z,q,r;F), I(z, q, r) = g(z, q, r; I). (69) 

The generating functions F and / satisfy the following 
system of equations, 



kff ~ 1 and F* < 1 in any case, 3 the value (|75|) is 
bounded from above by v + 1, which is not a large num- 
ber. Hence, the mean total number of H-regions of type 
F is at most v + 1 times larger than the mean total num- 
ber of bubbles of type F. (Both numbers are finite in fi- 
nite multiverses, and the relationship persists in the limit 
n to t -> oo.) 

Analogous results are obtained for regions of type /. 
We simply need to replace F with I and Kpp by Ku in 
Eq. (75]). 

The "Boltzmann brains" are created at a very small 
rate k bb per //-region of type j, while ordinary observers 
are created at a much larger rate per reheated 3-volumc 
HJ 5 in bubbles of the same type. We conclude that the 
RV-regulated abundance of "Boltzmann brains" is always 
negligible compared with the abundance of ordinary ob- 
servers in bubbles of the same type. 



(70) 
(71) 



F" — ZKpA + ZKpjI + 1'KppF, 

I~" = zkjb + zrqKjpF + Kill. 



These equations differ from the analogous Eqs. (|9"j)- ([lT 
in that the generating parameter q appears when a new 
bubble of type F is created, while the parameter r appears 
every time a new //-region of type F is created, which 
can happen via Hubble expansion of old F-bubbles as 
well as through nucleation of new F-bubbles. 

We would like to compare the mean number of H- 
regions of type F with the mean number of new bubbles 
of type F, so we compute (e.g. starting with /-bubbles) 



'n;I 



d"d r I 



(np 



d^dql 



(72) 



=0,r=l,q=l 



and take the limit n — -> oo. Taking the limit as n — > oo 
of the ratio d"F yr /d"F, q will yield the same result, but 
the limit of Eq. (f72")) is more straightforwardly analyzed. 

As before, we first obtain the equations for the deriva- 
tives F >q , I tq , F >r , I, r at q = r = 1 from Eqs. ([70]) - ([7T|) . 
The derivatives F iq and I iq satisfy the same equations as 
before, namely Eq. (|2"D)) , while F r and I. r satisfy 



M(z) 



' F. r ' 




kpfF 


T ,r _ 




ZKlpF 



whose solution is 



F, 



M-\z) 



kpfF 
zkifF 



(73) 



(74) 



Using the same arguments as in the previous section, we 
evaluate the limit of Eq. (172")) , 



dl l d r I 



n-^oo d^d q I 



= VF* " KFF + 1 • 



(75) 



IV. A GENERAL LANDSCAPE 

In the previous section we performed computations in 
a simple toy model of the landscape. Let us now consider 
a general landscape containing N vacua, labeled j = 1, 
N . The dimensionless transition rates Kj_ between 
vacua j and k are considered known. The main task is to 
compute the RV-regulated ratio of abundances of bubbles 
of kinds j and k. We will also compare the abundances 
of ordinary observers with that of Boltzmann brains. 



A. Bubble abundances 

It is convenient to denote by T the set of terminal 
bubble types and to relabel the vacua such that j = 1, 
N r are the "recyclable" (nonterminal) vacua. Thus, 
T = {N r + 1, N}. We start by considering the proba- 
bility p(n, {rij} ; k) of having nj bubbles of type j, with 
total n = rij bubbles of all types, to the future of 
an initial bubble of type k ^ T. The generating func- 
tion g(z, {qj} ; k) for this probability distribution can be 
defined by a straightforward generalization of Eq. , 



](z,{qi};k) 



E 

n,{m}>0 



p(n,{ ni }; k)z n Htf. (76) 



Here z is the generating parameter for n, and qj are the 
generating parameters for rij. The generating function 
can be written symbolically as the average 



g(z,{qi};k)^ {z n qT-q n N N 



n<oo;fc ' 



(77) 



where the subscript (n < oo) indicates that only the 
events with a finite total number of bubbles contribute 
to the statistical average. 



z=0,r= l,q— 1 



To analyze this simple result, we do not actually need to 
use the complicated decision procedure of Sec. IIII1 Since 



3 Typically F» < 1 but we can do with a weaker bound F, < 1 
here. 
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For terminal bubble types k £ T, the definition (175)) 
yields g(z, {qj} ; k) = 1 since there are no further bubbles 
to the future of terminal bubbles, thus n = rij 
probability 1. 

The generating function g(z, {qj} ;k), k = 1, . 
isfies the system of N r nonlinear equations 



with 



, N r sat- 



g»(z,{qj};k) 



E 

+ E Zq '' 



zqiKk->ig(z,{qj};i) 



Kk^i + Kk^kg{z, {qj} ; k), 

(78) 



where v = e 3 as before, and the quantities n k ^ k defined 
by Eq. {6]). A derivation of Eqs. (|78| will be given below 
in Sec. ITVDl 

The solution of Eqs. (|78|) can be visualized as N r an- 



alytic functions g(z, {qj} ; k), k = l,...,N r , of the free 
parameters z and qj (j = 1,...,N). Arguments similar 
to those of Sec. Mil show that the physically significant 
solution of Eqs. (|75|) is the "main branch" that has the 
following asymptotic form at z — ► 0, 



g(z;k) 



+ 0(z 



1v-\\ 



(79) 



(see also the argument at the end of Sec. IIV Dp . 

Once the generating function g(z, {qj} ; k) is deter- 
mined, the RV measure gives the ratio of the mean num- 
ber of bubbles of type j to that of type k as 



lim 



p(k) n^oo d™d qk g(z, {qi};i') 



(80) 



z=0,{ 9i }=l 



where, for clarity, we wrote explicitly the type %' of the 
initial bubble. We will now use the methods developed in 
Sec. IIII Al to reduce Eq. (|80| to an expression that does 
not contain limits and so can be analyzed more easily. It 
will then become evident that the limit (|50|) is indepen- 
dent of i'. 

For a fixed bubble type j, we first consider the deriva- 
tives d qj g(z, {qt} ;k),k = 1, N r , evaluated at {qi} = 1. 
Let us denote these N r derivatives by hj(z; k), 

h 3 (z;k)=^- g(z,{q t };k), k=l,...,N r . (81) 

These quantities are conveniently represented by an N r - 
dimensional vector, which we will denote by \hj{z)) using 
the Dirac notation (although no connection to quantum 
mechanics is present here). This vector satisfies an inho- 
mogeneous linear equation that follows straightforwardly 
by taking the derivative d qj at {qi} = 1 of Eqs. ([78]) . That 
equation can be written in the matrix form as follows, 



JV 

E 



M ki (z)hj(z;i) = zKk^jg(z; j), k ^ j, (82) 



JV 

E 



M ki {z)g, qj {.z; i) = 0, k = j, 



(83) 



where we denoted by g(z;j) = g(z;{qi} = the solu- 
tion of Eqs. ([75]) at {<&} = 1, written as 

g»(z;k) = z ^2 K k^ig(z;i) 

i£TA=ik 

+ z^K k ^i + K k ^ k g(z;k), (84) 

iGT 

while the matrix M k i(z) is defined by 
M{z)^M kl {z) = \^~ 1{z]k) - Kk ^ k -~ l ) (85) 

\-ZK k ->i, kf=t. 



For convenience we rewrite Eqs. (|82[) (|83|) in a more con- 
cise form, 



M(z)\h j (z)) = \Q j (z)), 



(86) 



where \Qj) is the vector with the components \Qj) i , i 
1, N r given by the right-hand sides of Eqs 



\Qj(z)) t = zn t ^g{z; j) [1 - Sij] . (87) 
The solution of Eq. (|86[) is found symbolically as 

M- 1 (z)\Q J (z)), 



\hj(z)) 



provided that the inverse matrix Af _1 (z) exists. Of 
course, it is impractical to obtain the inverse matrix ex- 
plicitly; but we will never need to do that. 

The next task is to determine z for which the matrix 
M(z) remains nondegenerate, detM(z) ^ 0. It will be 
shown in Sec. IIV El that there exists an eigenvalue Xo(z) 
of M(z) such that Ao(z*) = at some real z» > 0, while 
Xo(z) > for z < z*. Moreover, the (left and right) eigen- 
vectors corresponding to the eigenvalue \o(z) are always 
nondegenerate and can be chosen with all positive com- 
ponents. At the same time, all the other eigenvalues of 
M(z) remain nonzero for all z < z*. It will also be shown 
(see Sec. IIV Fl) that the following asymptotic expansion 
holds near z = z*, 



Xq(z) = ciy/z* — z + 0(z te — z), C\ > 0. 



(89) 



One can see from the small-z asymptotics (|79|) and from 
Eq. (|55|) that, for small enough z, M(z) contains very 
large positive numbers on the diagonal and very small 
negative numbers off the diagonal. Hence detM(z) > 
for all sufficiently small z. Since the eigenvalues of M(z) 
remain nonzero except for \o(z) that first vanishes at 
z = z„, we conclude that detAf(^) remains positive for 
all < z < z*, and that z — z* is the smallest positive 
value of z for which det M(z) = 0. 

Now we restrict our attention to < z < z* , for which 
detM(z) > and M^ 1 (z) exists. We need to analyze 
the behavior of M _1 near z w z*. We treat M(z) as 
an operator in A^-dimensional real space V. The matrix 
M is not symmetric and may not be diagonalizable. In- 
stead of diagonalizing M(z), we split the space V into 
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the 1-dimensional subspace corresponding to the small- 
est eigenvalue \o{z), and into the (N r — l)-dimensional 
complement subspace V\. We know that the eigenvalue 
Ao(z) is nondcgenerate. Thus, we may symbolically write 



M(z) = A (z) \v°(z)) (u°(z)\ + Mi(z), 



(90) 



where and (u°| are the right and the left eigenvectors 
corresponding to Xq(z), normalized such that 



(u°\v°) = 1, 



(91) 



and it is implied that M\{z) vanishes on and (u°\ 
but is nonsingular on the complement space V\ . In other 
words, we have 

M \v°) = X \v ) , (u°\M=Xq(u°\, (92) 
Mi\v°) = 0, (u°|Mi = 0. (93) 

There exists a matrix M~,h < that acts as the inverse to 

Mi when restricted to the subspace V\ and again vanishes 
on and (u°\. Using that matrix, we may write the 

inverse matrix M _1 explicitly as 



M~ 



X ^ = TTT + M i(y t )( z )- 

Ao(z) ^ u 

Hence, the solution (|88|) can be written as 



(94) 



\hj(z)) 



X (z) 



Now it is clear from Eq. 

z = z* as oc (z* — z) 



-1/2 



\v°)+M-f Vi) \Q j ). (95) 

that \hj(z)) diverges at 
It also follows that \hj(z)) 



does not diverge at any smaller real z. 

We can also show that Xq(z) cannot vanish at some 
complex value of z that is closer to z = than z = z*. If 
Xo(zl) = with a complex- valued then a derivative of 
the generating function, such as d z g(z;j), would diverge 
at z = z'^. Using the definition (|76p of g and substituting 
qi = 1, we find that the following sum diverges, 



dg{z;j) 



dz 



E 

n>l 



(4)" 1 np(n;j) 



(96) 



The sum of absolute values is not smaller than the above, 
and hence also diverges: 



dg(z; j) 



dz 



z —\ z * I 



(97) 



n>l 



So d z g has also a singularity at a real value z = \z^\. As 
we have shown, z — z* is the smallest such real-valued 
singularity point; hence \z'*\ > z*. It follows that z* 
is equal to the radius of convergence of the series (|96|) . 
which is a Taylor series for the function d z g. There re- 
mains the possibility that a singularity z^ is located di- 
rectly on the circle of convergence, so that \z'J\ = z* and 



zi — z^e 1 ^ with < (f> < 2ir. This possibility can be ex- 
cluded using the following argument. The function d z g 
can have only finitely many singularities on the circle 
| z | = z* ; infinitely many singularities on the circle would 
indicate an accumulation point which would be an essen- 
tial singularity, i.e. not a branch point. However, by con- 
struction g{z; j) is an algebraic function of z that cannot 
have singularities other than branch points. Since g(z;j) 
is regular at z = z*, all branch points of g must be of 
the form (z — z^) s with s > 0. Using the explicit for- 
mula ([23)) for the n-th derivative of an analytic function 
with a branch cut singularity of the form (z — z^) , we 
find that the coefficients np(n;j) of the Taylor series ([96]) 
decay at large n asymptotically as 

np(n;j) = -d n z d z g ex F( " 4 ' ]_~ " ] «p (l + 0{n^)) 
n! 1 (n + 1) 



ocn- s [zi)~ n (l + Oin- 1 )) 



(98) 



The function g has a finite number of singularities z^ = 
z^e 1 ^ on the circle \z\ — z*, and each singularity gives a 
contribution of the form (|98p . Hence, we may estimate 
(for sufficiently large n, say for n > no) 

np(n;j) = c rT s z: n (l + O^- 1 )) , c > 0. (99) 



Then the partial sum for n > no of the series (|96j) is 
estimated by 



n>no 



E 

n>no 



c n 



~ E ! 

n>no 



(100) 



The latter series converges for s > and < <j> < 2tt. 
(The neglected terms of order n _1_s build an absolutely 
convergent series and hence introduce an arbitrarily small 
error into the estimate.) Therefore, the series (|96|) also 
converges at z* ^ z* , contradicting the assumption that 
other singularities exist on the circle |z| = z*. 

We conclude from these arguments that the singularity 
of g(z;j) nearest to z — in the complex z plane is 
indeed at a real value z = z„, and all other singular 
points z^, satisfy the strict inequality \z'J\ > z». The same 
statement about the locations of singularities holds for 
the generating functions g(z, {q^} ;j) and hence for their 
derivatives such as \hj{z)). 

To compute the final expression ([50)1 . we need to eval- 
uate the n-th derivative d z of the functions hj(z;i) at 
z = and for very large n. To this end, we use the for- 
mula (|25[) . which requires to know the location z = z* of 
the singularity of hj(z;i) nearest to z = 0. We have just 
found that this singularity is at a real value z = z* and 
has the form 



hj(z;i) 



1 



■(u°\Qj)vl z^z r , (101) 



Cly/zT 

where v® is the i-th component of the eigenvector |i; ). 
The value of the proportionality constant c\ is not re- 
quired for computing the ratios ([501) . Using Eqs. ([2"S")) . 
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([57]) . and (jlOip . we evaluate the limit as 



B. Example landscape 



p(k) (u°\Q k )v° 



Hijtj^Z* 



(102) 

Since the component u2 cancels, we find that the 
limit ([50)) is indeed independent of the initial bubble type 
i'. We also note that the normalization of the eigenvector 
(u°\ is irrelevant for the ratio. 



The formula (|102[) can be simplified further for 
nonterminal types j, k if we use the relationship 
(u a (z*)\M{z*) = together with the explicit defini- 
tion (153 of M: 



= -^2ui(z*)z*Ki-+j + u° J (z„) ^g"^ 1 (z^;j) - Kj^j 

(103) 



The last term in the square brackets in Eq. (|103|) can 
be neglected since Kj-,j ~ 1 while g(z*;j) <C 1. Hence, 
Eq. (TlU2"]) is simplified to 



PC?) _ 



g( z *;j) 



l/u 



(104) 



This is the main formula for the RV-regulated relative 
abundances of bubbles of arbitrary (nonterminal) types 
j and fc. For a terminal type k £ T, one needs to use 
Eq. (fTU2"|) together with ,g(z; k) = 1. 

We note that the expression (1104[) depends on the com- 
ponents u° of the left eigenvector of the matrix M(z*) 
with eigenvalue Ao(z*) = 0. Although we have been 
able to evaluate the limit (|80|) analytically and obtained 
Eq. (|104j) , the task of computing the values of 2* , g(z* ; j) , 
and w?(z*) remains quite difficult. The outline of the re- 
quired computations is as follows: One first needs to de- 
termine g(z; j) (j = 1, N r ) as the "main branch" of the 
solution of Eqs. (154")) with the small- z asymptotic given by 
Eq. (|79p . The functions g(z;j) determine the matrix ele- 
ments Mij (z) using Eq. (jSSJ) . Then one needs to find the 
smallest value z = z* > such that the determinant of 
the matrix M(z) vanishes. Finally, one needs to compute 
a left eigenvector (u°(^*)| that corresponds to the eigen- 
value Ao(z*) = of the matrix M(z*). The mathematical 
construction shown below guarantees that z* exists and 
that (u°(z*)| is nondegenerate and has all positive com- 
ponents; the results are independent of the normalization 



of (t 



However, a brute-force numerical computation 



of these quantities appears to be impossible due to the 
huge number N r of the simultaneous equations (|84[) and 
to the wide range of numerical values of the coefficients 

> j m a typical landscape. Even the numerical value 

of z* is likely to be too large to be represented efficiently 
in computers. In the next section we will consider an 
example landscape where an analytic approximation can 
be found. 



We begin with some qualitative considerations regard- 
ing the behavior of the functions g(z; k). 

To determine g(z*;k), we need to follow the main 
branch g(z; k) as the value of z is increased from z = 
until det M(z) vanishes, which will determine the value 
z = z*. We note that the asymptotic form (|75)l is valid 
in some range near z = 0. Since g{z; k) <C 1 at those 
z, the matrix M(z) is dominated by large positive diag- 
onal terms g~~ 1 (z;k)5 k j. As z increases to z„, all the 
functions g(z;j) also increase while detM(z) decreases 
monotonically to zero (this statement is proved rigor- 
ously in Sec. HVE|) . One can visualize the changes in 
the matrix elements of M(z) if one notes that the pos- 
itive diagonal terms decrease with z while the negative 



off-diagonal terms (Mj k 



*fej 3 ^ k) grow in mag- 



nitude. Eventually det M(z) vanishes at z — z* such that 
the off-diagonal terms become sufficiently large at least 
in some rows and columns of the matrix M(z). 

Let us determine an upper bound on z such that 
Eq. (|79p remains a good approximation for the gener- 
ating function g(z;j) = g(z,{qi — This will be 
the case if the term z ^ i&T in the right hand side 
of Eq. (|84|) dominates over all other terms, for every k. 
Denoting for brevity 



(105) 



the total transition probability from k to all terminal 
vacua, we have therefore the conditions (for every k) 

ZK k ^ T > g(z;k), ZK k _> T > z ^ K k _ig(z;i). 

(106) 

These conditions are satisfied, consistently with Eq. (fT9"|) . 
if z is bounded (for every k) simultaneously by 



Z <C K fe ^T' z ^ 



E 

i^T,i^k 



K k^ri v 

<- 

Kk^T 



(107) 



Since typically the rate of transitions to terminal vacua 
is much smaller than the rate of transitions to dS vacua, 
one can expect that this regime will include z = 1. Nev- 
ertheless, one expects that det M(z) remains nonzero for 
these z. The detailed consideration of the FABI model 
in Sec. IIII Al showed that the value of z* lies outside the 
regime (|107[) unless the landscape parameters are fine- 
tuned. In a general landscape, one expects the analo- 
gous fine-tuning to be much stronger or even impossible 
to satisfy. In other words, one expects that the func- 
tions g(z*; k) at least for some k will violate the condi- 
tions (|106p . although these conditions might be satisfied 
for other k. 

Now we would like to estimate the value of z for which 
the matrix M(z) first becomes degenerate. We can use a 
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general theorem for estimating the eigenvalues of a ma- 
trix (Gershgorin's theorem, see (59[, chapter 7). For each 
of the diagonal elements Mkk, k — l,...,N r one needs 
to draw a circle in the complex A plane, centered at the 
diagonal element Mkk with radius 



Pk 



Ei M *i- 

i^tk 



Thus one obtains N r circles 



|A- M kk \ < Pk, 



1. 



(108) 



(109) 



called the Gershgorin circles. The Gershgorin theorem 
says that all the eigenvalues of a matrix Mjk are located 
within the union of these circles. An elementary proof is 
as follows. An eigenvector Vi with eigenvalue A can be 
normalized such that its component of largest absolute 
value is equal to 1. Let v\ = 1 be this component (renum- 
bering the components if necessary), then it is easy to 
show that A is within a circle with center M%x and radius 
pi . Namely, we start from the eigenvalue equation, 



E 



(110) 



and use the properties v\ = 1 and \vi\ < 1 for i ^ 1: 



lA-MnI 



MuVi 



<Ei m «i 



Pi- 



(111) 



Applying the theorem to the matrix M(z) at small z such 
that Eq. (|79| is a good approximation for the generating 
function g(z;j), we find that the Gershgorin circles are 
centered at large positive values 



-gu 



-\z-k) 



K-k—>k 



1 

V {zK k - 



» 1, 



while the radii of the circles are small, 



Pk= Z E K k- 



< 1. 



(112) 



(113) 



Hence, for small z none of the circles can contain A = 0, 
so the matrix M{z) is nondegenerate. As z increases, 
the radii pk increase while the centers Mkk move mono- 
tonically towards zero. The circle closest to A = is 
the one with the center closest to zero and the largest 
radius. This circle is labeled by k with largest rates 
Kk—>T and Kk—n ■ This value of k corresponds to the high- 
energy vacua in the landscape, for which tunneling to 
any other vacua is much easier than for low-energy vacua. 
Therefore the eigenvalue \{z») — will be located inside 
the Gershgorin circle(s) centered at Mkk{z*) for k corre- 
sponding to high-energy vacua. The Gershgorin circle to 
which the eigenvalue belongs indicates the largest com- 
ponent of the eigenvector. Hence, one expects that the 



eigenvector (u°(z*)| has the largest values of its compo- 
nents u° corresponding to the high-energy vacua k. The 
Gershgorin theorem requires that Ao(z*) = be inside 
the circle centered at Mkk(z*) with radius pk{z*), i.e. 



M k k(z*) = -g* 1 (z*;k) - Kk^k < z* E 

i=£k,i<£T 



Kfc- 



(114) 

Disregarding the terms Kk—>k, which are small in com- 
parison with the remaining terms, and using Eq. (|84p for 
<?(z*; k) we obtain the following condition, 



zZ > 



(115) 



Proceeding further requires at least an estimate of g(z*] i) 
for all i, which is so far not available. 

While we are as yet unable to obtain an explicit es- 
timate of z* and for an arbitrary landscape, let us 
consider a workable example that is more realistic than 
the FABI model. In this example there is a single high- 
energy vacuum, labeled k = 1, and a large number of 
dS low-energy vacua, k — 2, ...,N r , as well as a number 
of terminal vacua, k = N r + 1,...,N. We assume that 
the downward tunneling rates Ki^i for i ^ 1 are much 
larger than the rates Ki^i or Ki—,j for i,j ^= 1. Then 
g(z; 1) ^> g(z;i) for i 7^ 1, and so the first Gershgorin 
circle is the one closest to A = 0. Hence we may expect 
that the eigenvalue \o(z) always belongs to that circle. 
Moreover, it is likely that z* is within the regime p06[) 
for the low-energy vacua (k ^ 1) but not for the high- 
energy vacuum k = 1. We will now proceed with the 
calculation and later check that these assumptions are 
self-consistent. 

If the regime (|106[) holds for every low-energy vacuum 
k =/= 1, we have 



j(z;k) w z v K v h _ 



k = 2,...,N r . 



For the high-energy vacuum k — 1 we use Eq. 
rectly to find 



(116) 
(H di- 



g(z; i) = 



Ki^ T + z E *) + K!->ig(z; 1) 

E K l^i K 



+v 



(117) 



where we have neglected the terms with k,\^t arL d «i_>i. 
The condition that the regime (| 1 06[) holds for k ^ 1 is 



+ t > z E K 'i 



jg {z;j), i = 2,...,N r . (118) 



In the matrix M(z) the off-diagonal elements Mjj for 
=/= 1, j =/= 1 are much smaller than all other elements. 
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Hence, we can approximate M{z) by a matrix that has 
only its first row, first column, and the diagonal elements, 



M(z) 



( M n 



V 



zkn 



-3«l->2 






-ZKl^3 





-ZKl->N r \ 
' 



M 



N r N r 



M kk » -gi 1 {z-k). 



(119) 



This approximate matrix is much easier to analyze; in 
particular, its determinant and the eigenvectors can be 
computed in closed form. The determinant of this matrix 
is 



detM(z) w Mn...M NrNr 



Z Ki^iKi^i 



Therefore, the condition det M(z* 



§ Mn(z)M u (z) 

(120) 

can be written as 



«1- 



tKi-fl 



i=2 



M;;(z*) 



(121) 



Using Eqs. fTTSjl . (fTT7| . and fTl"9"|) . we transform this 
condition to 



' N r 

E 

,i=2 



EI/-1 

j=2 



(122) 

However, one does not need this value apart from check- 
ing explicitly that the assumptions (|118p hold. 

The left eigenvector ^u (z*)| corresponding to the 

eigenvalue of the matrix M(z*) is found approximately 
as 



1. 



Z*Ki^,k 



(123) 



where the normalization u\ = 1 was chosen arbitrarily 
for convenience. A perturbative improvement of this ap- 
proximation along the lines of Refs. 0, H3, IH, |4t| may 
be possible, but the precision obtained from the present 
approximation is sufficient for our purposes. 

The bubble abundance ratios (|104p for nonterminal 
types are then expressed as follows, 



pU) 

p(k) 

P (k) 



, V-lJV r v 

1 Z^i=2 K l^i K i^T 



j,k = 2,...,N r , 

k = 2,...,N r 



(124) 
(125) 



These equations are the main result of this section, yield- 
ing RV-regulated bubble abundances for a landscape with 
a large number of low-energy vacua. The formula Q125P 
agrees with that obtained in Sec. IIII Al for the ratio 
p(F) : p(I) 1/l> in the FABI landscape, which may 
be considered a special case of the present model with 

N r = 2. 



C. "Boltzmann brains" 

We now investigate the abundance of "Boltzmann 
brains" in a general landscape, relative to the abundance 
of ordinary observers. 

We first need to derive the equations for the suitable 
generating functions, analogously to Eqs. (|70p - (|7Tj) . Let 
us introduce the generating function 



9{z,{qi} ,{n};j) 



, (126) 



n tot <°°;j 



where rij is the total number of bubbles of type i and Ni 
is the total number of if-regions of type i. The number 
of BBs in bubbles of type i is proportional to Ni with 
a proportionality constant nf B , which is the (extremely 
small) nucleation rate of a BB per Hubble 4- volume H~ . 

The generating function g(z, {q{\ , {r{\ ; j), which we 
will denote for brevity by g(...,j), satisfies a system of 
equations analogous to Eq. (|75j) . 

9"(-,j) = 'Y^Kj^kzq k r k g{-,k) + K j -> j r j g(...,j). 

(127) 

Let us compute the RV-regulated ratio of the number of 
-ff-regions of type j to the number of bubbles of the same 
type j. This ratio is given by 



lim 



lim 



d%d rj g{...,i') 



(128) 



z=0,r i = l,g i = l 



where i' is the initial bubble type. To evaluate the 
limit, we use the methods developed in Sec. IIV Al The 
derivatives d q .g = \hj(z)), as defined in Eq. t[5Tj). satisfy 
Eq. (j5B"|) . The vector of derivatives 



_d_ 



g{...,k) = \pj(z))t 



(129) 



=1,9; 



satisfies the linear equations that follow from Eq. (|127[) . 

M(z)\p j (z)) = \(3 j (z)), (130) 

where \/3j) is the vector with the components \j3j) v i = 
1 , . . . , N r defined by 

\/3 j (z)) i = zK^g(z;j) [1 - <%] + Kj^jg(z;j)6ij. (131) 



By the same considerations that lead to Eq. (I102j) . we 
now obtain 



lim 



d2d rj g(...,i>) 



00 d™d qj g(...,i') 



<«°l&> 



z=0,ri=l,qi 



( u °\Qj 



(132) 

where, as before, is the unique eigenvector of the 

matrix M(z*) with eigenvalue 0. (The normalization of 
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m 



is again irrelevant.) Using the definitions ([57)1 and 
} and the relationship p03[) , we compute 



(u°\Qj 



v9~ 



(133) 



The last ratio is always very close to 1 since typically 
g(z*;j) <C 1 while Kj—>j ~ 1. In particular, g(z*;j) <C 1 
for vacua j with low-energy Hubble scale, since for those 
vacua we may approximate g(z*;j) by Eq. (f79|) . 



l(z*,j) 



■ ^2 K j~>i 



< 1. 



(134) 



Therefore, the RV-regulated ratio Nj/rij of the total 
number of ff-regions of type j to the total number of 
bubbles of type j is never large. 

Using this result, we can now estimate the RV- 
regulated ratio of BBs to ordinary observers. The number 
of ordinary observers per one .ff-region of type j is not 
precisely known but is presumably at least of order 1 or 
larger, as long as bubbles of type j are compatible with 
life. On the other hand, the number of Boltzmann brains 
per iJ-region, i.e. within a four-volume HJ A of space- 
time, is negligibly small. It follows that the abundance 
of BBs in the RV measure is always negligible relative to 
the abundance of ordinary observers in the same bubble 
type. 



D. Derivation of Eq. (TTgj) 

To derive Eq. (78), we need to consider the expan- 
sion of a single initial H -region during one Hubble time. 
Within the "inflation in a box" model, an initial H-region 
of type j is split after one Hubble time 5t = HJ 1 into 
v = e 3 statistically independent "daughter" f/-regions. 4 
Each of the daughter iJ-regions may change its vacuum 
type from j to k =/= j (k = l,...,N) with probability 
Kj^k- The quantity Kj—>j was defined for convenience by 
Eq. © to be the probability of not changing the bubble 
type j during one Hubble time. 

The generating function g is defined by Eq. (|76p. 

g(z,{ qi };j) = (z^qT- q n N N ) ntal<00 , 3 (135) 

where rii is the number of bubbles of type i, while the 
notation (... )„ tot<oc .j stands for a probabilistic average 



4 To avoid considering a non-integer number v of daughter regions, 
we may temporarily assume that v is an integer parameter. At 



the end of the derivation, we will set v 
equations. 
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evaluated for the initial iJ-region of type j on the sub- 
ensemble of finite total number of bubbles n to t ■ Note that 
the initial bubble of type j is not counted in n, or n to t- 
Our goal is to obtain a relationship between g(z, {qt} ;j) 
and the generating functions g(z, {q{\ ; k) with k ^ j. 

To this end, we equate two expressions for the av- 
erage (z" tot T7, ;<??*) ^ ■■ The first expression is the 
left-hand side of Eq. (|135|1 . The second expression is 
found by considering the v daughter i7-regions evolved 
out of the initial f/-region and by using the fact that the 
same average for a daughter region of type k is equal to 
g{z, {q{\ ; k). However, two details need to be accounted 
for: First, the generating functions g(z, {<?;} ; k) evalu- 
ated for the daughter regions do not count the daughter 
bubbles themselves. Second, the type k of each of the 
daughter regions is a random quantity. Before deriving a 
general relationship, let us illustrate the procedure using 
an example. 

It is possible that, say, only two of the daughter regions 
change their type to k while all other daughter regions 
retain the initial bubble type j. Denote temporarily by 
Pkkj...j the probability of this event; binomial combina- 
torics yields 

Pkkj...j = 2\(y- 2j\ K ^f° K j-?j- ( 136 ) 
Then the average of z ntat q^ receives a contribution 
Pkkj...jZ 2 ql [g(z,{ qi };k)f [g(z, {q l } ; j)] 2 ^ (137) 

from this event. The factor z 2 describes two additional 
bubbles that contribute to rttot; the factor q\ accounts 
for two additional bubbles of type k; no factors of qj 
appear since no additional bubbles of type j are gener- 
ated. Finally, the powers of g account for all the bubbles 
generated in the daughter £f-regions, but these generat- 
ing functions do not count the daughter 77-regions them- 
selves. Those daughter 77-regions are explicitly counted 
by the extra factors z 2 q\. 

To compute the average, we need to add the contribu- 
tions from all the possible events of this kind. Since all of 
the daughter ii-regions are independent and statistically 



equivalent, the average (£ nta * FT. q^] 



"tot<oo;j 



splits into 



the product of v averages, each evaluated over a single 
daughter region. 

The average over a single daughter region of type j has 
contribution from transitions to other types k ^ j and 
a contribution from the event of no transition. Let us 
first consider a terminal type k. With probability Kj^k 
a given daughter region becomes a vacuum of type k. 
Thereafter, no more bubbles will be nucleated inside it; 
the average of (z ntat Yl-q^*) „ - over that daughter 
region is simply zqk, meaning that there is a total of 
one bubbles and only one bubble of type k. Hence, the 
contribution of that event to the statistical average is 
Kj^kzqk' 

Now let us consider a nonterminal type k =/= j, 
k (I T. The corresponding contribution to the av- 
erage is Kj^kzqkg(z, {qi} ; k). Since we have defined 
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g(z, {qi} ; k)= 1 when k is a terminal bubble type, we 
may write the contribution as Kj^zq^giz, {qi} ; k) for 
both terminal and nonterminal types k =/= j. 

Finally, we consider the case of k = j (the daughter 
region retains the original bubble type). Since no new 
bubbles were nucleated, the contribution to the average 
is simply Kj^jg(z, {qi} ;j) without any factors of z or q k . 

Putting these ingredients together, we obtain an equa- 
tion for g(z,{qi};j), 

g( z ,{li}'J) 

^2 Kj^kzq k g(z, {qi} ; k) + Kj^jg(z, {q t } ; j) 

(138) 

This is equivalent to Eq. (|78|). One can also verify that 
the binomial expansion of Eq. (|138[) indeed yields all the 
terms such as the one given in Eq. (|13T[) . 

We can now analyze the behavior of g(z, {qi}; j) in 
the limit z — > 0. By definition, the generating function 
g(z, {qi} ; j) is a power series in z whose coefficient at z n 
is equal to the probability of the event that a multiverse 
has exactly n bubbles to the future of the initial bub- 
ble j. It is clear that this power series starts with the 
term z v , corresponding to the probability that the initial 
bubble j expands exactly into v terminal bubbles, sig- 
nalling the global end of the multiverse. The probability 
of having fewer than v bubbles in the entire multiverse 
is equal to zero. 5 The next term of the binomial expan- 
sion is of order z 2ly_1 since it is the product of z v ~ Y and 
g itself. Therefore, the small-z behavior of the gener- 
ating function g(z, {qi} ;j) must be given by Eq. ([75)) . 
This condition, together with analyticity in z, selects the 
unique physically relevant solution of Eqs. f|T8|) . 



E. Eigenvalues of M(z) 



According to the definition (|8"5|) . the matrix M{z) has 
positive elements on the diagonal and nonpositive ele- 
ments off the diagonal. Such a matrix can be rewritten 
in the form 



M(z) = fil - A(z), 



(139) 



5 The property that there are exactly v daughter bubbles is, of 
course, an artifact of the "inflation in a box" approximation. In 
the actual multiverse, one has bubbles of spherical shape that can 
intersect in complicated ways, so a given bubble may end in one, 
two, or any other number of terminal bubbles. However, we are 
using the box approximation to obtain results in the limit of very 
large total number of bubbles, so we disregard the imprecision 
in the description of multiverses with a very small total number 
of bubbles. Also, the value of v may be considered a variable 
parameter of the "box" model; the final results will not be overly 
sensitive to the value of v. 



where a constant fi > is introduced, the notation 1 
stands for an identity matrix, and A(z) is a suitable non- 
negative matrix, i.e. a matrix with all nonnegative ele- 
ments. For instance, we may choose /i as the largest of 
the diagonal elements of M. The theory of nonnegative 
matrices gives powerful results for the eigenvalues of ma- 
trices such as M and A (see e.g. the book [5!|). For the 
present case, the most important are the properties of 
the algebraically smallest eigenvalue of the matrix M. 

It will be convenient to drop temporarily the argument 
z since all the results of matrix theory will hold for every 
fixed z. By the Perron- Frobenius theorem (see chap- 
ter 9) , under the condition of irreducibility 6 a nonnega- 
tive matrix A has a unique nondegenerate, real eigenvalue 
ccq > such that all the eigenvalues of A (which may be 
complex- valued) are located within the circle |A| < ao 
in the complex A plane. This "dominant" eigenvalue ao 
has a corresponding (right) eigenvector |t>°) that can be 
chosen with all strictly positive components v® > 0. The 
same property holds for the relevant left eigenvector (u°\ 
(the matrix need not be symmetric, so the left and the 
right eigenvectors do not, in general, coincide). Therefore 
it is possible to choose the eigenvectors (ti°| and \v°) such 
that the normalization (u°\v°) = 1 holds. This normal- 
ization will be convenient for further calculations, and so 
we assume that such eigenvectors have been chosen. 

It follows that \v°) and are also the right and left 
eigenvectors of the matrix M with the eigenvalue 



\q = fi — a , 



(140) 



while all the other eigenvalues of M are located within 
the circle — A| < ao in the complex A plane. Since 
ao > 0, all the other eigenvalues of M are strictly to 
the right (in the complex plane) of the real eigenvalue 
Ao- In other words, Ao is the eigenvalue of M with the 
algebraically smallest real part. 

Restoring now the argument z of the matrix M, we 
find that M always has a real, nondegenerate eigenvalue 
Xo{z), which is at the same time the eigenvalue with the 
algebraically smallest real part among all the eigenval- 
ues of M(z). We know that detM(z) > for suffi- 
ciently small z; hence Xo(z) > for those z. Moreover, 
detM(z) will remain positive as long as An (2) > 0, since 
no other eigenvalue can become negative unless Ao (z) first 
becomes negative. We will now show that det M(z) can- 
not remain positive for all real z > 0. It will then follow 



6 The irreducibility condition means that any two recyclable vacua 
in the landscape can be connected by a chain of transitions with 
nonzero nucleation rates. This condition has been discussed in 
Refs. |23t [Sill . If some subset of vacua form a "disconnected 
island" in the landscape, such that transitions to and from the 
"island" are forbidden, one can regard the "island" as a separate 
irreducible landscape and apply the same technique to it. Hence, 
we consider only irreducible landscapes in this work. 
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by continuity of Xq(z) that there will be a value z» such 
that X (z) > for all < z < z* but Aq(z*) = 0. 

We will use the property that the inverse matrix 
Af _1 (z) has all positive elements as long as Xq(z) > 
(equivalently if ao < /x). The derivation of this property 
is simple: 



and hence 



M" 



{pi -A) = fj,- 1 l + fi-^A + f j,- 6 A J + ..., (141) 



which yields explicitly a matrix with all nonnegative ele- 
ments. [The matrix- valued series in Eq. (|141[) converges 
because all the eigenvalues of A are strictly smaller than 
fi by absolute value.] Moreover, the irreducibility con- 
dition means that some chain of transitions will connect 
every pair of recyclable vacua; this is equivalent to saying 
that for any vacua i, j there exists some integer s such 
that A s has a nonzero matrix element (A 8 )^. Hence, ev- 
ery matrix element of M~ l is strictly positive as long as 
AoCz) >0. 

Further, we can deduce that g(z;j) is a strictly in- 
creasing, real-valued function of z for those z for which 
Xq(z) > 0. To show this, we consider the vector \d z g) 
whose components are the N r derivatives d z g(z;j), j 
1 . .... Y, . It follows from Eq. ((751) that the vector \d z g) 
satisfies the inhomogeneous equation 

M(z) \d z g) = |0 , (142) 
where we denoted by |£) the vector with the components 
Ck{z) =^n k ^ig{z;i). (143) 

i^k 



The solution of Eq. (fUg) is 



\d z g)=M-\z)Uz)) 



(144) 



Since all the matrix elements of M _1 (z) are positive 
and all the components of |£) are nonnegative as long 
as g(z; i) > 0, it follows that all the components of \d z g) 
are strictly positive. Equation (|75)l shows that g(z: i) > 
for sufficiently small z > 0, and it follows that g(z; i) will 
remain positive for all z > such that Xq(z) > 0. There- 
fore, g(z; i), i = 1, N r are strictly increasing functions 
of z for all those z. 

Nevertheless, the functions g(z;i) are bounded from 
above. To see this, consider the relationship 



(u°\M = A (? 



(145) 



written in components as 



u%z) 



1 A-l, .x 

-g v 



- X! u i( z ) ZK i^j = Ao^°. 

(146) 

Since all the components m° are strictly positive (as long 
as Xq(z) > 0), it follows that 



(147) 



g(z;j) < {vKj^j) - 1 < 



1 1 



(148) 



Furthermore, we can show that Xq{z) monotonically 
decreases as z grows. This follows from the perturbation 
theory formula for nondegenerate eigenvalues, which al- 
lows us to express dX^/dz as a matrix product with nor- 
malized left and right eigenvectors, 

^l = {u o {z) \^L\ v o {z) y (149) 

As we have just shown, g(z;i) grows with growing z, so 
dM/dz is a matrix with all nonpositive elements. Since 
the vectors and have strictly positive compo- 
nents while at least some matrix elements of the nonpos- 
itive matrix dM/dz are strictly negative, we obtain the 
strict inequality 

dXn(z) i n , . I dM i n , .> 

-^ 1 = (u°(z)\ — \v°(z))<0. (150) 

Similarly, we can show that det M(z) monotonically de- 
creases with z: 



d detM(z) = (detM(z))Tr(M- 1 ^) < 



dz 



dz 



(151) 



since it was already found that the matrix M 1 (z) has 
all positive elements while dM/dz has all nonpositive el- 
ements. However, the monotonic decrease alone of Ao(.z) 
and of det M(z) is not yet sufficient to establish that the 
matrix M(z) actually becomes singular at some finite z. 

The results derived so far — the monotonic behavior 
of g(z; j) and Xq(z), the positivity of the matrix elements 
of M _1 , the bounds on j — hold for all z for which 
Xo(z) > 0. Now we will show that Ao(z) cannot remain 
positive for all real z > 0. We can rewrite Eq. (|84|) as 



zVit^,j(z;i) = g»(z;k) - n k ^ k g(z;k). (152) 



Using the property g(z;j) > and the bound (|148|) . we 
obtain (for every k) an upper bound on z, 



g-(z;k) - K k ^ k g{z;k) 



E 



i^T.i^k 



n k ^ig(z;i) + n k - 



(153) 

In other words, no real-valued solutions of Eq. (|84| ex- 
ist for larger z. Let us then show that an upper bound 
on z contradicts the assumption that Xo(z) > for all 
z. We know that there exists a real- valued solution 
branch g(z;j) near z = such that < g(z;j) < oo 
and < d z g(z;j) < oo for all those z > for which 
this solution branch remains real- valued. Hence, g(z;j) 
can be viewed as a solution of a differential equation 
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dzg(z~,j) = \d z g) with continuous coefficients and every- 
where positive right-hand side. The solution of such dif- 
ferential equations, if bounded, will exist for all z > 0. 
Indeed, if the solution g{z;j) existed only up to some 
z = z\, we would have, by assumption, Xq{z\) > and 
hence a finite value g(zi;j) > and a finite derivative 
d z g(z\;j) > 0. So the solution g(z;j) could then be con- 
tinued further to some z > Z\ . Therefore, the real- valued 
solution branch g(z;j) must exist for all z > 0. This is 
incompatible with the bound (|153D . 

We conclude that there exists a value z» > such that 
Ao(z*) = but Xq(z) > for all < z < z*. Within the 
range < z < z* the functions g(z;j) grow monotoni- 
cally but remain bounded by Eq. (I148|) . while Xq(z) and 
det M(z) both decrease monotonically to zero. 



F. The root of X (z) 

It remains to establish that Xq(z) indeed has the 
form (|89[) near z = z*. We again restrict our attention 
to the interval < z < z* where Xq(z) > 0. For these z 
we expand Xq(z) in Taylor series and express the value 
Xq(z^) = as 

= A (z») = Xo(z) + (z* -z) + 0[(z* - zf], 



hence 



dz 



dX (z) Xq{z) 



+ 0(z* - z). 



(154) 



(155) 



dz z* — z 

We then use Eq. (|149[) to express dXo /dz in another way, 



^ = (u°(z)\±M(z)\v°(z)) 



■.(u°{z)\ 



Y,^rr L -A9(z;i) + d z M(z) 

, dg(z;i) 



\v°(z)), 
(156) 



where in the second line we interpreted M(z) as a func- 
tion of N r variables g(z;i), i — l,...,N r , and explicitly 
of 2, in order to express d/dz through d/dg and d/dz. 
Using Eq. (|144j) . we then find 



Ao(*) 



Since we are only interested in the qualitative behavior 
of Ao(z) at z = z„, we do not need to keep track of the 
complicated coefficients in Eq. (jl5T|) . Near z — z* we 
have, by Eq. ([M]). 



1 



« )(«°|+O(l), 



(158) 



so the dominant singular terms in Eq. (|15T|) near z — z* 
are 



Xo(z) 



Ci 
X (z) 



0(1). 



(159) 



Therefore we obtain 



X (z) — ciVz* — z + 0(z* — z). 



(160) 



The positivity of Xo(z) for z < z* entails c\ > 0. This 
concludes the derivation of Eq. 
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